Spin Orientation and Spin Precession in Inversion- Asymmetric Quasi 
Two-Dimensional Electron Systems 



R. Winkler 

Institut fur Technische Physik III, Universitat Erlangen-Niirnberg, Staudtstr. 7, D-91058 Erlangen, Germany 

(Dated: May 14, 2003) 

Inversion asymmetry induced spin splitting of the electron states in quasi two-dimensional (2D) 
systems can be attributed to an effective magnetic field B which varies in magnitude and orientation 
as a function of the in-plane wave vector kn . Using a realistic 8x8 Kane model that fully takes into 
account spin splitting because of both bulk inversion asymmetry and structure inversion asymmetry 
we investigate the spin orientation and the effective field B for different configurations of a quasi 
2D electron system. It is shown that these quantities depend sensitively on the crystallographic 
direction in which the quasi 2D system was grown as well as on the magnitude and orientation of 
the in-plane wave vector ky . These results are used to discuss how spin-polarized electrons can 
precess in the field B(kii). As a specific example we consider Gao.47lno.53As-InP quantum wells. 



I. INTRODUCTION 

Spin degeneracy in a two-dimensional (2D) system is 
due to the combined effect of spatial inversion symme- 
try and time inversion symmetryi If the spatial inver- 
sion symmetry is lifted spin-orbit interaction gives rise 
to a spin splitting of the electron states even at a mag- 
netic field B = 0. In quasi 2D systems the B = spin 
splitting can be caused by the bulk inversion asymme- 
try (BIA) of the underlying crystal structure- as well as 
by the structure inversion asymmetry (SIA) due to, e.g., 
an electric field £ perpendicular to the plane of the 2D 
system. - The B — spin splitting is of considerable inter- 
est both because of its importance for our understanding 
of the fundamental properties of quasi 2D systems 4-5 - 6 - 7 
as well as because of possible applications in the field of 
spintronicsi^ 

Common III-V and II- VI semiconductors such as 
GaAs, InSb, and HgCdTe, have a zinc blende structure. 
To lowest order in the wave vector k BIA spin splitting 
in these systems is characterized by the so-called Dressel- 
haus termSi whereas spin splitting due to SIA is charac- 
terized by the Rashba term£ Often the discussion of spin 
splitting is restricted to these lowest-order terms ^ULILH 
Spin splitting of higher orders in k can be fully taken 
into account by the 8x8 Kane models or the 14 x 14 
extended Kane model^ The higher-order terms can be 
quite important for a quantitative discussion of B = 
spin splitting. iSiiii 

For a given in-plane wave vector k|| we can always 
find a spin axis (S(k||)) local in k|| space such that we 
have spin-up and spin-down eigenstates with respect to 
the axis (S(k||)). Note that we cannot call the spin- 
split branches i£±(k||) of the energy surface spin-up or 
spin-down because the direction of (S) varies as a func- 
tion of k|| such that averaged over all occupied states 
the branches contain equal contributions of up and down 
spinor components. This reflects the fact that in nonmag- 
netic materials we have at B = a vanishing magnetic 
moment. 

The spin orientation (S(k||)) can be attributed to an 



effective mag netic field B(ky) fRefs. HTl . A discussion 
of (S(k||)) based on the lowest-order terms in the effective 
spin-orbit interaction has previously been given by sev- 
eral authors, see, e.g., Refs. ll8ll9l20l2ll22l In the present 
paper we compare these results with our calculations of 
(S(k||)) and the field S(kii) using the more realistic 8x8 
Kane modeli^ that takes into account both SIA and BIA 
up to all orders in kii . It will be shown that for larger ky 
the higher-order terms result in important modifications 
of (S(k||)) and S(k||). 

Datta and Das have proposed a novel spin transistor— 
where the current modulation arises from the precession 
of spin-polarized electrons in the effective field S(kii), 
while ferromagnetic contacts are used to preferentially 
inject and detect specific spin orientations. Recently, ex- 
tensive research aiming at the realization of such a device 
has been under way^ Here we will use our results for 
the field 8(kii ) in order to discuss spin precession and its 
tunability for different device configurations. It will be 
shown that for certain configurations the precession of 
spin-polarized electrons is determined only by the tun- 
able SIA spin splitting; but it is essentially independent 
of the magnitude of BIA spin splitting. For other config- 
urations the tunability of spin precession is significantly 
suppressed due to the interplay of SIA and BIA. 

We would like to emphasize that the present results 
apply only to electrons with an (effective) spin j = 1/2. 
Holes in the topmost valence band, on the other hand, 
have an effective spin j = 3/2 fRef. l24l) . Therefore, spin 
orientation and spin precession in quasi 2D hole systems 
is qualitatively different from spin orientation and spin 
precession in quasi 2D electron systems. Hole systems 
will thus be covered in a future publication. 



II. SPIN ORIENTATION OF 2D ELECTRON 
STATES 

In the following we want to discuss the wave vector 
dependent spin orientation (S(k||)) for different models of 
spin splitting. We will compare the analytical results for 



the Rashba model and Dresselhaus model with our more 
realistic calculations based on the 8x8 Kane Hamiltonian 
that takes into account SIA and BIA spin splitting up to 



all orders in k 



General Discussion 



First we want to discuss the spin orientation in the 
presence of SIA. Here to lowest order in the in-plane wave 
vector k|| = (k x ,k y ,0) the spin splitting is characterized 
by the Rashba Hamiltonian 3 



-ffsiA = a (a x k y - <jyk x ) , 



(1) 



where a x and a y are Pauli spin matrices and a is a pref- 
actor that depends on the constituting materials and 
on the geometry of the quasi 2D system. If we use 
polar coordinates for the in-plane wave vector, k|| = 
fcii (cos ip, sin ip, 0), the spin splitting is given by 



£| IA (k||) 



±afc 



independent of the angle <p and the eigenstates are 



l^ 1A (k||)> 
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with r 1 1 = (x,y,0) and envelope functions £k|, (2). In 
Eq. J2J we have assumed that the Rashba coefficient a 
is positive. The spin orientation of the eigenstates © is 
given by the expectation value with respect to the vector 
<j of Pauli spin matrices. 



(<r(k||)> ee <V±(k||)|<r|V±(k||)) 



(4a) 




Note that Eq. J3J is independent of the envelope function 
£k|i (z) and the magnitude k\\ of the in-plane wave vector. 
The spin orientation (J3J of the eigenfunctions J2J as a 
function of the direction of the in-plane wave vector is 
indicated by arrows in Fig. [Ja) • 

Next we want to discuss the spin orientation in the 
presence of BIA spin splitting. For quasi 2D systems 
in a quantum well (QW) grown in the crystallographic 
direction [001] the Dresselhaus term becomesi2*ii 

H BIA = 77 [a x k x (k 2 y - (k 2 z )) + a y k y ((k 2 z ) - *£)] (5) 

with a material-specific coefficient r\. This equation can 
easily be diagonalized. We obtain a spin splitting 



£f A (k||) = ±rjk\J{kl 
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FIG. 1: Lowest order spin orientation (<r) of the eigenstates 
|^>±(k||)) in the presence of (a) SIA and (b) BIA. The inner 
(outer) circle shows (tr) along contours of constant energy for 
the upper (lower) branch E+ (E-) of the spin-split dispersion. 



We see here that in leading order of k|| the Dresselhaus 
term gives rise to a spin splitting independent of 
the direction of k|| that is apparently very similar to 
the Rashba spin splitting J2J. Nevertheless, the corre- 
sponding wave functions are qualitatively different due 
to the different symmetries of the terms {Q and (J5J. If 
we neglect the terms cubic in k|| the eigenfunctions in 
the presence of Dresselhaus spin splitting are 



l^ 1A (k||)) = 



2tt 



1 



so that 



M k n)>. 




(7) 



(8) 



The spin orientation © of the eigenfunctions as a 
function of the direction of the in-plane wave vector is 
indicated by arrows in Fig. ^b). For the Rashba spin 
splitting we see in Fig. ^a) that if we are moving clock- 
wise on a contour of constant energy S(kii) the spin vec- 
tor is rotating in the same direction, consistent with the 
axial symmetry of the Rashba term. On the other hand, 
Eq. (jHJ and Fig. Cfb) show that in the presence of BIA 
the spin vector is rotating counterclockwise for a clock- 
wise motion in k|| space. 

In the above discussion we have assumed that the wave 
functions are two-component spinors. In general, the 
quasi 2D eigenstates of a multiband Hamiltonian are of 
the form2£ 



.ikiirii 



Wkii)) = X^« k n( z ) u "( r ) 



(9) 



with envelope functions £ n k||(-z), and u n (r) denotes the 
band edge Bloch function of the nth bulk band. Here we 
must evaluate the expectation value of 



i 



orb i 



(10) 
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where the identity operator l or b refers to the orbital part 
of |^(k||)). For the 8x8 Kane modeli^ containing the 
bands T§, fg, and we obtain for i = x,y,z 



Si = 



2 7 

3 1 

-2U] 
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(11) 



where Ji denotes the matrices for angular momentum 
j = 3/2, and the matrices Ui are defined in Ref.^J Once 
again the expectation value (^|S|^) is a three-component 
vector that can be identified with the spin orientation 
of the multicomponent wave function \ip). We remark 
that while the vector (a) of a spin 1/2 system is always 
strictly normalized to unity, this condition is in general 
not fulfilled for the spin expectation value (S) of multi- 
component single particle states. This is due to the fact 
that in the presence of spin-orbit interaction we cannot 
factorize the multicomponent wave function @ into an 
orbital part and a spin part. However, for electrons the 
deviation of |(S)| from unity is rather small (typically 
less than 1%) so that it is neglected here. 

For free electrons in the presence of an external mag- 
netic field B the unit vector (er) is parallel to the vector 
B. Following this picture we can attribute the B — spin 
splitting in quasi 2D systems to an effective magnetic field 
B(k||) parallel to (S(k||)). Obviously the magnitude of 
this effective magnetic field should be related to the mag- 
nitude of the B — spin splitting. However, depending 
on the particular problem of interest it can be convenient 
to define the magnitude of spin splitting in two different 
ways: The energy difference AE = E + (k^) — £L(ky) 
characterizes the magnitude of spin splitting for a given 
wave vector k|| whereas the wave vector difference Ak 
characterizes the magnitude of spin splitting at a fixed 
energy E. While the former is relevant, e.g., for Ra- 
man experiments^ the latter quantity is an important 
parameter, e.g., for spin relaxation&ii and for the spin 
transistor proposed by Datta and Dasii& 

In the following we want to explore the second def- 
inition where the effective magnetic field is given by 
B = (S)Afc. Our precise definition of Ak is illustrated 
in Fig. |21 For the given energy E and a fixed direction 
<p of the in-plane wave vector k|| = feii(cos</?, sin</?,0) we 
determine k|| =F Ak/2 such that E = £+(ky - Ak/2) = 
E-(k\\ +Ak/2). Here E + {EJ) denotes the upper (lower) 
branch of the spin-split dispersion. Then we define 



B = (S)+ Afc = -(S)_ Afc 



(12) 



with the sign convention that the field B is parallel to the 
effective field felt by the electrons in the upper branch 
J?+(kii) and we have used the short-hand notation 

(S)± - (V±(k|| T Ak/2) | S | V±(k|| T Ak/2)). (13) 

We remark that for a parabolic band with effective mass 
m* plus Rashba term (JTJ the wave vector difference Afc 




+Ak/2 



FIG. 2: For the given energy E and a fixed direction of the 
in-plane wave vector ky we determine kii =p Ak/2 such that 
E = E+(k\\ - Ak/2) = £7_(k|| + Ak/2). Here E+ (E-) 
denotes the upper (lower) branch of the spin-split dispersion. 



can be evaluated analytically*^ 



Ak 



Rashba 



2m* a 



(14) 



independent of the magnitude of fey . From an exper- 
imental point of view it should be kept in mind that 
spin splitting is often measured by an alyzing Sh ubnikov- 
de Haas oscillations, see, e.g., Refs. I2fiu27ll28ll2fl Such 
experiments yield spin subband densities N± which are 
directly related to Ak 



Ak = 



(15) 



provided we can ignore anisotropic contributions to 
B = spin splitting. (However, see also Refs. 

mm) 

The definition l|12|) presupposes that the spin expecta- 
tion values (S)+ and (S)_ are strictly antiparallel to each 
other. In Eq. Q we saw that for the Rashba Hamiltonian 
this condition is fulfilled exactly. This is closely related 
to the fact that for the Rashba Hamiltonian the spin sub- 
band eigenstates |i/rp A (k||)) and |^L IA (ky)) are orthogo- 
nal - independent of the magnitude of ky and ky as long 
as the wave vectors ky and ky are parallel to each other ml 
In general, |V>+(ky - Ak/2)) and |V-(ky + Ak/2)) are 
only approximately orthogonal so that (S)+ and (S)_ 
are only approximately antiparallel. However, we find 
that the angle between the vectors (S)+ and (S)_ is al- 
ways very close to 180° with an error < 1° so that we 
neglect this point in the remaining discussion. 

Even though we can evaluate the spin expectation 
value (S) for each spin subband separately we do not 
attempt to define an effective magnetic field B for each 
spin subband. This is due to the fact that B is commonly 
used to discuss phenomena like spin relaxation&i£ and 
spin precessioniS (see below) which cannot be analyzed 
for each spin subband individually. 

The allowed directions of the effective magnetic field B 
can readily be deduced from the symmetry of the QW. 
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The spin-split states for a fixed wave vector ky are or- 
thogonal to each other, i.e., the spin vectors of these 
states are antiparallel. The spin orientation of eigen- 
states for different wave vectors in the star of k|| arc 
connected by the symmetry operations of the system^ 
Accordingly, only those spin orientations of the spin-split 
eigenstates are permissible for which every symmetry op- 
eration maps orthogonal states onto orthogonal states. In 
a QW grown in the crystallographic direction [001] the 
effective field B is parallel to the plane of the quasi 2D 
system. Indeed, the field 13 due to SIA is always in the 
plane of the well. For growth directions other than [001], 
the effective field due to BIA has, however, also an out-of- 
plane component. In particular, a symmetric QW grown 
in the crystallographic direction [110] has the point group 
Ci v w^ Here the BIA induced field S(ky) must be perpen- 
dicular to the plane of the QW (to all orders in ky). This 
situation is remarkable because D'yakonov-Perel' spin re- 
laxation is suppressed if the spins are oriented perpendic- 
ular to the 2D plane^& Note also that in [110] grown 
QWs 3 vanishes for k|| || [001] because here the group of 
k|| is C'2v which has merely one irreducible double group 
representation, Ts, which is two-dimensionaL -- 

B. Numerical Results 

The analytically solvable models Q and (JSJ allow one 
to study the qualitative trends of BIA and SIA spin split- 
ting in quasi 2D systems. The largest spin splitting can 
be achieved in narrow-gap semiconductors where the sub- 
band dispersion is highly nonparabolic. Therefore, we 
present next numerically calculated results for B(ky) ob- 
tained by means of an accurate 8x8 Kane Hamiltonian 
(Tg, Tg, and Tj) including off-diagonal remote band con- 
tributions of second order in k (Refs. ll3Ul6T) . First we 
analyze BIA spin splitting that is always present in zinc 
blende QWs. In Fig. Ota) we show the effective field 
(|12fl along contours of constant energy for a symmetric 
GaAs QW grown in the crystallographic direction [001] 
with a well width of 100 A. The dimensions of the arrows 
in Fig. |3 are proportional to \B\ = Ak. We remark that 
typical Fermi wave vectors of quasi 2D systems are of the 
order of the in-plane wave vectors covered in Fig. [3J 

For small in-plane wave vectors k|| the effective field in 
Fig. Eta) is well described by Eq. JSJ. For larger wave vec- 
tors the effective field becomes strongly dependent on the 



direction of k||. In particular, we see that for ky || [110] 
the effective field reverses its direction when we increase 
. This reversal reflects the breakdown of the linear ap- 
proximation in Eq. (JSJ). For wider wells this breakdown 
occurs at even smaller wave vectors fey, consistent with 
Eq. ©. 

More specifically, Eq. © predicts for k|| || [110] a re- 
versal of the direction of B(k||) when k 2 = 2(fc 2 ), in- 
dependent of the material specific coefficient r\. Note, 
however, that (k 2 ) depends on the material specific band 
offset at the interfaces. For the system in Fig. Of a) we 
find in good agreement with Eq. © that the reversal of 
B(k||) occurs for k\\ « y / 2(fcf) » 0.029 A -1 . For compar- 
ison, we show in Fig. Etb) the effective field S(ky) for a 
symmetric Gao.47Ino.53As QW with the same well width 
100 A like in Fig. Eta). Even though BIA spin splitting- 
is smaller in Gao.47Ino.53As than in GaAs, higher-order 
corrections are more important in Gao.47Ino.53As due to 
the smaller fundamental gap of this material. Here we 
have (kf) « 3.6 x 10~ 4 A" 2 so that v / 2(k 2 ) w 0.027 A" 1 . 
On the other hand, the reversal of the direction of B(ky) 
occurs for fey w 0.021 A -1 . This illustrates the effect of 
higher orders in BIA spin splitting that were neglected 
in Eq. (J5J) but fully taken into account in the numerical 
calculations in Fig. El [Note that in Fig. EI a ) the effective 
field 3 has been amplified by a factor of 50 whereas in 
Fig. Elb) it has been amplified by a factor of 100.] 

Gao.47Ino.53As QWs can have a significant Rashba 
spin splitting^ so that these systems are of interest for re- 
alizing the spin transistor proposed by Datta and Das>i£ 
In Fig. Etc) we show the effective field B(kii) for the 
same well like in Fig. Etb) assuming that we have SIA 
spin splitting due to an electric field £ — 20 kV/cm, 
but all tetrahedral terms that give rise to BIA spin split- 
ting were neglected. The numerical results are in good 
agreement with what one expects according to Eqs. (J1J 
and ifHjl. Figure Etd) shows the effective field S(kii) 
for a Gao.47Ino.53As QW when we have both BIA and 
SIA spin splitting. Due to the vectorial character of B 
we have regions in ky space where the contributions of 
BIA and SIA are additive whereas in other regions the 
spin splitting decreases due to the interplay of BIA and 
SIA. This is consistent with the well-known fact that in 
the presence of both BIA and SIA the spin splitting is 
anisotropic even in the lowest order of k|| (Ref. Using 
Eqs. Q and JSJ we obtain 



^BIA+SIA 



n o? + a V (k 2 - 2(fc 2 )) sim» + r, 2 (fc 2 ) 2 + (±fc 2 - (k 2 z ))k 2 sin(2^) 



± fc|, yta 2 - 2ar](k 2 z ) sin(2^) + r] 2 (k 2 ) 2 ± O(kjj). 

I 



(16a) 
(16b) 



In Figs.|3ta-d) we have considered QWs grown in the crystallographic direction [001] so that the effective field 
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GaAs [001], BIA, £ = 0kV/cm Ga.47ln.53As [001], BIA, £ = 0kV/cm Ga.47ln.53As [001], no BIA, £ = 20 kV/cm 




k x (A" 1 ) k x (X" 1 ) k x (V) 

Ga.47ln.53As [001], BIA, £ = 20kV/cm 



T 1 1 1 r 




k x (A" 1 ) 

FIG. 3: Effective magnetic field S(kii) for (a) a GaAs-Alo.3Gao.7As QW and (b-e) a Gao.47lno.53As-InP QW, both with a well 
width of 100 A. In (a), (b), and (e) we assume that we have a symmetric well with BIA spin splitting only, (c) shows B(kii) due 
to an external field of £ = 20 kV/cm but neglecting BIA while (d) shows jB(kri) when we have both BIA and SIA spin splitting 
(again for £ — 20 kV/cm). While (a-d) refers to a QW grown in the crystallographic direction [001] we have assumed in (e) 
that the QW was grown in [110] direction. The dimensions of the arrows are proportional to \B\ = Ak. For Gao.47Ino.53As, we 
have amplified S(kii) by a factor of 100, for GaAs it has been scaled by a factor of 50. All calculations are bas ed on an 8x8 
Kane Hamiltonian (Fg, Tg, and T^) including off-diagonal remote band contributions of second order in k fRef. llJTtl) . 



S(kii) is always in the plane of the QW. For comparison, 
we show in Fig. Ete) the effective field B(k||) for a sym- 
metric Gao.47Ino.53As QW grown in the crystallographic 
direction [110] with k x || [00T] and k y || [TlO]. Here S(k||) 
is perpendicular to the plane of the QW. 3 - For asymmet- 
ric QW's grown in the crystallographic direction [110] the 
effective field S(kjj) is given by a superposition of an in- 
plane field as in Fig. 02c) and a perpendicular field as in 
Fig. He). 



III. SPIN PRECESSION OF 2D ELECTRON 
STATES 

A. Datta Spin Transistor 

We want to briefly recapitulate the mode of opera- 
tion of the spin transistor proposed by Datta and Dasi& 
(see Fig. 0J. We assume that the semiconducting chan- 
nel between the ferromagnetic contacts is pointing in 
x direction, i.e., electrons travel with a wave vector 
k|| = (k x ,0,0) from source to drain. A gate in z direc- 
tion gives rise to a tunable Rashba coefficient a. In this 
subsection we want to ignore the Dresselhaus spin split- 
ting JSJ. When the spin-polarized electrons in the ferro- 
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FIG. 4: Qualitative sketch of a Datta spin transistor,' 18 Black 
arrows indicate the spin polarization in the ferromagnetic con- 
tacts (FM) and the semiconducting channel (white). Gray 
arrows indicate the effective magnetic field B{k x ) in the semi- 
conducting channel. A top gate is used to tune the spin pre- 
cession by applying an electric field £ perpendicular to the 
semiconducting channel. 



magnetic source contact are injected at x = into the 
semiconducting channel we must expand its wave func- 
tion \ipi) in terms of the spin-split eigenstates \ip± A (k x )}. 
Here it is the basic idea of the spin transistor that the po- 
larization of the electrons in the source contact is chosen 
perpendicular to B{k x ) — (0,6^,0). The states \tpi) thus 
contain equal contributions of the spin-split eigenstates 
\ip± A (k x )) . Assuming that the electrons in the source 
contact are polarized in +z direction we get (neglecting 
the envelope functions ^ky (z) which are unimportant in 
the present discussion) 



= o)) = 



i 



(17a) 



The basis states on the right hand side of Eq. I|17a|) prop- 
agate with wave vectors k x =p Afc /2 as depicted in Fig. [21 



|-0iO)) = exp [i (fc x - Afc/2) x] 



+ exp[z (k x + Ak/2)x] 



(17b) 



Due to the different phase velocities of the basis states in 
Eq. 117bl) we thus get 



sin(— Afc x) 
(S(,))=| 

cos(Afc x) 



(18) 



This equation can be visualized by saying that the spin 
vector (S) of the state \ipi(x)) precesses around the effec- 
tive field B(k x ) = (0,B y ,0) (see Fig. gj. Note, however, 
that conventional spin precession^ takes place as a func- 
tion of time t whereas in Eq. (|18fl the spin precesses as a 
function of position x. 

If finally the drain contact at x = L is ferromagnetic, 
too, the electrons can exit the semiconducting channel 



only if the spin orientation (S(x = L)) of the electrons 
matches the polarization Prj of the drain contact, 



cosx 



(S(x = L)) , 



(19) 



where x denotes the angle between Pd and (S(x = L)). 
A large positive value of cosx indicates that the elec- 
trons can easily exit the semiconducting channel whereas 
a large negative value indicates that the spin-polarized 
current is suppressed. Assuming that Ps,d I! [001] we 
obtain from Eq. ljl"%|l 



cosx = cos(Afc L). 



(20) 



A tunable device is achieved if the wave vector difference 
Afc is varied by changing the Rashba coefficient a, see 
Eq. dJ} and Fig. H 

In the above qualitative discussion we have ig- 
nored details such as the resistance mismatch at the 
interfaces^S^i^Si which are important for the practical 
realization of such a device. But these aspects do not af- 
fect the spin precession inside the semiconducting chan- 
nel which is the subject of the present investigation. 



B. Precession in the presence of BIA and SIA 

In the preceeding subsection we have assumed that 
only the Rashba term Q contributes to spin splitting. 
Here the effective magnetic field S(kii) that character- 
izes the spin orientation of the eigenstates © is always 
perpendicular to the direction k|| of propagation in the 
spin transistor. In general, we have both SIA and BIA 
spin splitting so that the effective field S(kii) is a more 
complicated function of k|| , see Fig- EI An arbitrarily ori- 
ented effective field B(k||) can be characterized by polar 
angles 9 and 4>, i.e., B = Afc [sin 9 cos </>, sin 9 sin <j>, cos 9] . 
The corresponding orthonormal eigenstates are 



IT) = 
II) = 



e"^/ 2 cos(#/2) \ 
e ^/ 2 sin(0/2) J 

_ e -#/ 2 sin (0/ 2 ) 
e 4 ^ 2 cos(0/2) 



(21a) 
(21b) 



For any values of the angles 9 and 0, the spin states 
(|21ll represent a basis of the spin 1/2 space. Similar to 
Eq. I|17a|) we can thus expand the wave function \ipi) of 
the spin-polarized electrons in the ferromagnetic source 
contact in terms of the basis states (|21|) 



\ipi(x = 0)) = cosu|T) 



'ID 



(22a) 



with angles u and v. Thus we get for the precessing 
electrons inside the channel 



\ipi(x)) = exp [i (k x — Afc/2) x] cos it |t) 

+ exp[i(fc a . + Afc/2)x] sin Me™ ||). 



(22b) 
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Then the overlap of the spin vector (S(x)) with the field 
& is given by 



3 • (SO)) = Afc cos(2u) 



(23) 



independent of the position x inside the channel and in- 
dependent of the phase e m . This equation shows that 
in generalization of Eq. (|18|) the spin is precessing on 
a cone around the effective field 3 where the cone an- 
gle is 2u. The precession amplitude Afc cos(2u) is the 
largest when u = 7r/4 so that in Eq. (|22J) we have equal 
contributions of the spin-split states ||) and |J,). This 
corresponds to the situation that the spin polarization 
Ps = (S(x = 0)) in the ferromagnetic source contact is 
perpendicular to S(kii). Spin precession is suppressed 
for u = and u — tt/2 when the spin polarization Ps in 
the ferromagnetic source contact is parallel to B(ki|) so 
that only one spin state l|21|l contributes in Eq. I|22|) . 

We have seen in Fig. 01 that for a fixed wave vector 
ky the orientation of S(kii) can change when the Rashba 
spin-orbit interaction is tuned by means of an external 
gate. It follows that the basic operating principle of the 
Datta spin transistor remains valid for the more general 
eigenstates (|21fl provided the polarization Pg of the fer- 
romagnetic source contact is orthogonal to B(ky) for all 
values of the external "knob" that is used to tune the 
spin-orbit interaction. If the condition Ps J- B(kii) is 
not strictly fulfilled the tunability of the spin transistor 
is reduced. We note that these conclusions are valid also 
for the more general eigenstates (JjJJ. 



C. Numerical Results 

We present next numerically calculated results for the 
spin precession in a spin transistor obtained by means 
of an 8 x 8 Kane Hamiltoniani2ii& that takes fully into 
account both BIA and SIA. According to Fig. |5]the ef- 
fective fields S(kii) due to BIA and SIA in a [001]-grown 
QW are always parallel to each other for k|| || [110] and 
k|| || [110]. On the other hand, for k|| || [100] the fields are 
perpendicular to each other so that we want to focus on 
these two extreme cases. We will again consider a 100 A 
wide Gao.47lno.53As-InP QW, and we assume that the 
distance between source and drain contact is L = 5 /an. 
For ease of notation we will use a suitably rotated co- 
ordinate system (Fig. QJ such that the electrons always 
propagate in x direction, i.e., kii = (fc x ,0, 0). We assume 
that the Rashba spin-orbit coupling is tuned by applying 
an electric field £ perpendicular to the plane of the quasi 
2D system (Fig.gJ. 

In Fig. we show the overlap cos x between the spin 
vector (S(ai = L)) and the polarization Pd of the drain 
contact as a function of electric field £ . We consider 
different polarization states Ps of the source contact and 
it is assumed that Ps || Pd- The results in Fig. can 
readily be understood by means of Fig. [3J (i) If k|| |( 
[110] or kii || [110] BIA is of little importance because 
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FIG. 5: Overlap cosx between the spin vector (S(x = L)) 
and the polarization Pd of the drain contact as a function of 
electric field £ in a 100 A wide Gao.47ino.53AsTnP QW with 
a channel length of L — 5 /im. In (a) we assume Ps,d || 
ky whereas in (b) we assume Ps,d || [001]. Different line 
styles correspond to different crystallographic directions of 
ky as indicated. The calculations are based on an 8 x 8 Kane 
Hamiltonian (F§, Tg, and T?) including off- diagona l remote 
band contributions of second order in k (Ref. Il^llfil ). 
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FIG. 6: Spin splitting Afc as a function of electric field £ in 
a 100 A wide Gao.47ino.53As-InP QW. We consider different 
magnitudes and different crystallographic directions of ky as 
indicated in the figure. 



Sbia II SsiA- Consistent with Eq. (|2U|) we thus get a 
sinusoidal dependence of cos x on £ with the same angle 
X for P S , D || k|| and P s ,d || [001], see Figs. EJa) and 
(b). We note that for fixed magnitudes of k|| and £ the 
angle x can be adjusted by changing the length L of the 
channel. In the present work L has not been optimized. 
Note that the smaller is the length L the larger must 
be the modulation of £ for switching the device, (ii) 
For Ps,d || k|| || [100] and £ = the spin precession 
is suppressed because Ps,d || B. In this case we have 
cosx = 1 independent of the channel length L. For £ > 
the spin states start to precess. Here spin precession 
and cos x are more complicated functions of £ because £ 
changes both the magnitude and orientation of B. (iii) 
For a QW grown in the high-symmetry crystallographic 
direction [001] the overlap cos x is symmetric with respect 
to £ > and £ < 0. In the latter case the roles of 



110] and kii || [110] are reversed. 



It is interesting to compare Fig. [5] with the magni- 
tude of spin splitting Afc as a function of electric field £ 
(Fig. EJ • We see that Afc depends rather sensitively on 
both the magnitude and orientation of the wave vector 
k||. Nevertheless, we obtain in Fig. [5] the same modula- 
tion of the overlap cosx as a function of £ for k|| || [110] 
and kj| || [110], independent of the magnitude of k|| [apart 
from a constant phase shift xo(k||)]. This is due to the 
fact that the relevant quantity for the spin transistor is 
not the absolute value Afc of the spin splitting, but the 
variation d(Ak)/d£. We see in Fig. that the latter 
quantity depends much more weakly on the magnitude 
and orientation of k|| . Furthermore, it is advantageous 



that the orientation of S(k||) is independent of the mag- 
nitude of the tunable part of the spin-orbit interaction. 
It can be seen in Fig. [3] that this condition is fulfilled for 
k|| || [HO] and k,| || [110] but not for k,| || [100]. There- 
fore the modulation of cos x as a function of £ is more 
pronounced in the former case than for k|| || [100], even 
though in all cases the spin splitting Afc shows roughly 
the same field dependence d(Ak)/d£. 



D. Spin Precession and Spin Relaxation 

For the Datta spin transistor it is advantageous to have 
a small spin relaxation in the semiconducting channel 
because spin relaxation is competing with the controlled 
spin precession in the channel. Typically, the dominant 
mechanism for spin relaxation in 2D electron systems is 
the one proposed by D'yakonov and Perel' (DP)iiL 3 £ It 
can be viewed as a spin precession in the effective field B 
that is randomized because B changes when momentum 
scattering changes the wave vector k|| of the electrons. 
DP spin relaxation can therefore be suppressed if (apart 
from a sign of B) the orientation of B is independent of 
the wave vector k|| and the spins of the propagating elec- 
trons are oriented parallel to B. Such a situation can be 
realized in a symmetric QW grown in the crystallographic 
direction [110] where B is perpendicular to the plane of 
the QW; 35 i 36 see Fig. |^e) . Similarly, in a QW grown in 
the crystallographic direction [001] with |a| = \i]\ we have 
in first order of k|| that B | [110] (or B | [110] depending 
on the sign of a and rf)w^ In both cases spin relaxation 
is suppressed only for a particular value of the Rashba 
spin-orbit coupling (i.e., a particular value of the field £). 
For the spin transistor it is preferable to have a regime of 
electric fields £ with suppressed spin relaxation so that 
we can switch between cos x = 1 and cos x = — !• 

Recently, an alternative spin transistor has been 
proposed^ that is less sensitive to spin relaxation. It 
uses the fact that not only DP spin relaxation can be sup- 
pressed if (apart from a sign of B) the orientation of the 
effective field B is the same for all wave vectors k|| ; but 
obviously spin precession is then suppressed, too. There- 
fore, if B || Ps,d electrons travel unperturbed through 
the device which corresponds to the "on" state. In a de- 
tuned system, on the other hand, B varies as a function of 
ky which implies that, in general, B |f Ps,d- Therefore, 
spin precession and/or DP spin relaxation reorient the 
spins in the channel. The spin vector (S(x = L)) thus no 
longer matches the polarization Pd of the drain contact 
so that the current through the device diminishes. 

Such a spin transistor can be built using a QW grown 
in the crystallographic direction [110]. Here it follows 
from Fig. I^e) that if the QW is symmetric then DP spin 
relaxation is suppressed because B is perpendicular to 
the plane of the QW for all in-plane wave vectors k|| . If 
Ps.d II B electrons thus travel unperturbed through the 
device. If the QW is made asymmetric by applying an 
electric field £ perpendicular to the plane of the well, the 
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current diminishes because of the onset of spin precession 
and/or DP spin relaxation. 

Alternatively, we can use a QW grown in the crystallo- 
graphic direction [001] (Ref . 201) . Here we can achieve in 
linear order of kn that 3 is independent of k| if \a\ = \rf\. 
This situation is approximately shown by the innermost 
contour in Fig. Old). Note, however, that higher orders 
in spin splitting [in particular the cubic term in Eq. <|16|l ] 
do not comply with the requirement that in the on state 
of the device the orientation of 3 should be independent 
of ky. Furthermore, we see in Fig. EUd) that only the 
orientation but not the magnitude of 3 is independent 
of k||. For electrons with ky || [110] we have actually 
B = whereas B is the largest for k|| |j [110]. In the 
former case (i.e., for Ps,d II k|| || [110]) changing £ di- 
minishes the current through the device because we have 
then Ps,d IBso that injected electrons precess around 
3. The electrons do not precess in the latter case be- 
cause we have Ps,d || 3 independent of £. DP spin re- 
laxation is highly anisotropic, too. Here the situation is 
actually reversed: We have large spin relaxation rates for 
those directi ons of (S) for which we have a large B(kii) 
(Refs. I43l44t) . Therefore, spin relaxation supports the 
switching of the device most effectively if k|| || [110]. We 
remark that an all-inclusive investigation of this question 
should explicitely evaluate spin relaxation as a function 
of k|| and (S). 

IV. CONCLUSIONS 

In general, the total B = spin splitting in inversion 
asymmetric 2D systems is determined by an interplay of 
spin splitting due to BIA, which is always present in sys- 
tems with a zinc blende structure, and the tunable spin 
splitting due to SIA. These spin splittings can be char- 
acterized by effective magnetic fields S(ky) that vary as 



a function of in- plane wave vector k||. The functional 
form of BsiA(k||) due to SIA is independent of the crys- 
tallographic direction in which a QW has been grown. 
Due to the axial symmetry of the Rashba term the field 
2?siA(k||) is always perpendicular to ky in the plane of 
the QW. Furthermore, it is only weakly dependent on the 
magnitude of ky. On the other hand, the field BBiA(k||) 
due to BIA depends sensitively both on the magnitude 
and orientation of ky as well as on the crystallographic 
direction in which the QW was grown. For QW's grown 
in the direction [001] the field iE?BiA(k||) is always in the 
plane of the QW whereas for QW's grown in the direc- 
tion [110] it is pointing perpendicular to the plain of the 
QW. For other growth directions the field SBiA(k||) has 
both in-plane and out-of-plane components. 

Electrons injected into a 2D semiconducting channel 
propagate with a certain in-plane wave vector k|| . If these 
electrons are spin-polarized such that the spinor \ipi) of 
the electrons is not a spin eigenstate of the system, the 
spin of the propagating electrons precesses in the effec- 
tive field B(k||). The precession is the largest if the spin 
orientation (S) of the electrons is perpendicular to the 
effective field B(kii). In a QW grown in the crystallo- 
graphic direction [001] it is thus advantageous that the 
electrons are injected in the in- plane directions [110] or 
[110] because here the field 2?(k||) is always perpendic- 
ular to the direction of propagation. For the direction 
[100], on the other hand, the fields due to BIA and SIA 
are perpendicular to each other so that the orientation of 
the total field S(k||) depends on the magnitude of BIA 
and SIA spin splitting. 
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